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Abstract. - Chaotic systems, that have a small Lyapunov exponent, do not obey the common 
random matrix theory predictions within a wide "weak quantum chaos" regime. This leads to a 
novel prediction for the rate of heating for cold atoms in optical billiards with vibrating walls. The 
Hamiltonian matrix of the driven system does not look like one from a Gaussian ensemble, but 
rather it is very sparse. This sparsity can be characterized by parameters s and g s that reflect the 
percentage of large elements, and their connectivity respectively. For g s we use a resistor network 
calculation that has direct relation to the semi-linear response characteristics of the system. 
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The heating of particles in a box with vibrating walls is 
a prototype problem for exploring the limitations of linear 
response theory (LRT) and the quantum-to-classical cor- 
respondence (QCC) principle. In the experimental arena 
this topic arises in the theory of nuclear friction [I] , and 
more recently in the studies of cold atoms that arc trapped 
in optical billiards [2J. It is also related to the analysis 
of mesoscopic conductance of ballistic rings [3] . Formally 
the dynamics is generated by a time dependent Hamilto- 
nian H[f(t)], where f(t) parametrizes the displacement of 
boundary, analogous to the time dependent electric field 
of the conductance problem. In typical circumstances the 
classical analysis predicts an absorption coefficient G that 
is determined by the Kubo formula [JHE]j leading to the 
"Wall formula" in the nuclear context, or to the analogous 
"Drudc formula" in the mesoscopic context. 

If upon quantization we get for the absorption coef- 
ficient an h dependent result, that does not correspond 
to the classical result, we call it an anomaly. The ques- 
tion arises what are the circumstances in which anomalies 
show up joT-112]. There are "microscopic circumstances" 
in which an anomaly is not a big surprise: (1) If f(t) 
is slowly varying, so-called quantum adiabatic parametric 
driving, then Landau-Zener transitions between neighbor- 
ing levels might be the dominant mechanism for heating 
[5], and hence QCC is not expected. (2) If f(t) is low 
frequency noisy driving, that induces Fcrmi-Goldcn-Rulc 
(FGR) transitions between neighboring levels only, the re- 
sult would be determined by the level spacing statistics, 
and hence QCC is not expected [TTj . 

In this Letter we identify a "weak quantum chaos 



regime" where a quantum anomaly shows up in quite typi- 
cal "mesoscopic circumstances" ', where QCC would be ex- 
pected by common-wisdom. 

Modeling.— Wc consider a weakly chaotic billiard that 
has linear size L and a convex wall of radius R. The 
Hamiltonian can be written schematically as 



H[f(t)] = n-f(t)F = Uo + U-f(t)F 



(1) 



Specifically with regard to the numerical example of Fig[T] 
TLo describes a non-deformed rectangular box of length 
L x = L = 1.5 (upper edge), and width L y = 1.0. The 
term U describes the deformation of the fixed (left) wall: 
it is an arc of radius R = 8 whose center of curvature is 
shifted upwards a vertical distance Ay = 0.1 to break the 
reflection symmetry. The term F is the perturbation due 
to the displacement f(t) of the moving (right) wall which 
can be regarded as a piston. Later we characterize the 
time dependence of /(£). 

Our interest is focused in circumstances in which the 
Lyapunov (correlation) time t R = R/v E is much longer 
than the ballistic time t L = L/v E , where v E = {2E/m) 1 / 2 
is the velocity of the particle. Turning to the quantum 
analysis we realize that the minimal model for % depends 
on two dimensionless parameters: 

u = L/ R [dimensionless deformation] (2) 

h = X E /L [dimensionless Planck const] (3) 

Here A E =27rftp lanck /(mz; E ) is the de Broglie wavelength. 
For a given deformation (R determines u) and energy win- 
dow (E determines K) we calculate the eigenvalues and 
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Fig. 1: Sketch of the billiard system of Eq. (fTJ). The unper- 
turbed billiard is a rectangle of size L x — 1.5 and L y — 1.0. The 
deformation U, due to the curvature of the left wall (radius 
R—8), is characterized by the parameter it = L y /R. In order 
to break the mirror symmetry the center of the curved wall is 
shifted upwards a vertical distance Ay = 0.1. The time depen- 
dent perturbation is due to the displacement f(t) of the right 
wall. In the numerics the units are chosen such that ftpi anck = l 
and the mass is m=l/2. The image in the background repre- 
sents the eigenstate -B n ~13618. 




Fig. 2: Image of the perturbation matrix. Image of 
the matrix X — {|F nm | 2 } for the billiard of Fig. 1 within the 
energy window 3500 < E n < 4000. This matrix is sparse. More 
generally it might have some texture. The latter term applies if 
the arrangement of the large elements is characterized by some 
pattern. 



cigenfunctions of H using the boundary clement method 
[13] . find the ordered eigenenergics E n , and calculate the 
matrix elements F nm using the formula 



f™ = ~ / \ {n \vW {m) {y)dy 

2m ./ 



(4) 



where tp( n ' (y) is the normal derivative of the nth eigen- 
function along the piston boundary. An image of a repre- 
sentative matrix is displayed in Fig|2j and its bandprofile 
is presented in FigJ3J 

The absorption coefficient.— Having in mind cold 
atoms in an optical trap, we regard the wall vibrations, 
say of the "piston" , as low frequency noisy driving. The 
power spectrum of f(t) is described by a spectral function 



S(u) = e 2 ^- 



6 ^— ex P 



(5) 



As is common in the mesoscopic context we assume its 
spectral support to be to c < l/t R , but larger compared 
with the mean level spacing. Accordingly, in the numerics 
it is natural to take uj c as matching the first minimum in 
the bandprofile of Fig £2 

Following [13] we assume that there are FGR tran- 
sitions between levels, whose rate is proportional to 
\F nm \ 2 S(E n —E m ). As a result the system absorbs energy 
in rate Gs 2 analogous to Joule heating. We define 
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This is the classical hard chaos result for the absorption 
coefficient, which is obtained, e.g. using a kinetic picture, 



if one neglects correlations between successive collisions. 
This is a straightforward adaptation of the well known 
"Wall formula" of nuclear physics, which is analogous to 
the "Drudc formula" in condensed matter physics. 

Objective.— Our objective is to calculate the actual 
absorption coefficient G, i.e. to go beyond the "Wall for- 
mula" prediction, taking into account the implications of 
having t R 3> t h , which is the case for small deformation 
(u <C 1). The calculation of the actual absorption coeffi- 
cient G will be done below either within the framework of 
LRT using the Kubo formula (getting G LRT ), or within the 
framework of semi-linear response theory (SLRT) [TTMT2] 
using a resistor- network calculation (getting G SL r T )- The 
correlations between collisions lead to an LRT result that 
we would like to write as G LRT = g c Go. Similarly it is 
convenient to write the outcome of the SLRT analysis as 
follows: 
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(7) 



If QCC considerations apply, then g s ~ 1 with small h 
dependent corrections. The LRT and SLRT numerical 
results for g arc displayed in Fig.^J and the details are 
presented in what follows. 

Conflicting expectations.— Both in LRT and in 
SLRT the result for G depends on the "average" over the 
near diagonal elements of \F nm \ 2 , i.e. those that are in 
the strip \E n — E m \ < uj c . The difference between LRT 
and SLRT is how this "average" is defined: as a simple 
algebraic average, or via a resistor network calculation. 
For a small deformation, first order perturbation theory 
(FOPT) implies that these couplings are oc u 2 . But as u 
becomes larger the common expectation, based on Wigner 
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Fig. 3: The band profile of the matrix, (a) The algebraic 
average and median along the diagonals of the X nm matrix 
versus ui = (E n —E m ). The vertical axis is normalized with re- 
spect to Coo, while the horizontal axis is u/v E - The classical 
power spectrum is presented to demonstrate the applicability 
of the semiclassical relation Eq. (J5|. The red line is the analyt- 
ical expression that applies to zero deformation. The quantum 
analysis is for R = 8 with 100 < E < 4000 (EW1), and with 
10000 < E < 14000 (EW2). The dotted vertical line is the fre- 
quency 1/£ L and the dashed one is l/f R . (b) Zoom of the 
u <C 1/£l region. For sake of comparison we display results 
also for R — 2. The vertical lines indicate the mean level spac- 
ing. The dashed red curves are a refined version of Eq. (|12[l . 



theory, is to have Lorcntzian mixing of levels, leading to 
oc 1/u 2 smearing. In the formally equivalent problem of 
a conductance calculation this implies G oc 1/u 2 , where 
u represents the strength of the disordered potential (in- 
stead of using the FGR or Wigner picture one can use the 
equivalent Drude picture where the Born mean free path 
is ex 1/u 2 ). On the other hand the semiclassical expecta- 
tion, based on kinetic consideration, is to have, because of 
the bouncing, enhanced energy absorption oc 1/u. Loosely 
speaking the latter expectation follows from the observa- 
tion that a sequence of 1/u correlated collisions with the 
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Fig. 4: SLRT vs LRT. The scaled absorption coefficient g c 
(LRT) and g — g s g c (SLRT) versus the dimensionless 1/h (up- 
per panel), and versus the dimensionless deformation parame- 
ter u = L/R (lower panel). Note that g = 1 is the prediction 
of the "Wall formula" , while the line is based on the classical 
analysis. In the upper panel the analysis has been done for the 
billiard of FigfT] The calculation of each point has been car- 
ried out on a 100 x 100 sub-matrix of X centered around the 
h implied energy E. The "untextured" data points are calcu- 
lated for an artificial random matrices with the same bandpro- 
file and sparsity (but no texture). The complementary lower 
panel is oriented to show the small u dependence. The anal- 
ysis is based on a truncated matrix representation of Ho + U, 
within an energy window that corresponds to 1/h ~ 9. Due 
to the truncation there is some quantitative inaccuracy with 
regard to the larger g values. 



piston is like a single big collision. The purpose of the fol- 
lowing paragraphs is to resolve this confusion by adopting 
a generalized random matrix theory (RMT) perspective. 

RMT modeling.— So called "quantum chaos" is the 
study of quantized chaotic systems. Assuming that the 
classical dynamics is fully chaotic, as in the case of a bil- 
liard with convex walls (Fig[T]), one expects the Hamilto- 
nian to be like a random matrix with elements that have 
a Gaussian distribution. This is of course a sloppy state- 
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mcnt, since any Hamiltonian is diagonal in some basis. 
The more precise statement is following |14j : Assume that 
H generates chaotic dynamics, and consider an observable 
F that has some classical correlation function C(t), with 
some correlation time i R . Then the matrix representation 
F nm in the basis of % looks like a random banded ma- 
trix. The bandwidth is h/t R . If t R is small, such that 
the bandwidth is large compared with the energy window 
of interest, then the matrix looks like it is taken from a 
Gaussian ensemble. 

What emerges in our numerical example, we would 
like to call "weak quantum chaos" (WQC) circumstances, 
for which the traditional RMT modeling does not apply. 
Namely, in such circumstances it is not enough to char- 
acterize F nm by its semiclassically-dctermined bandpro- 
file. Rather one should further characterize F nm by its 
quantum-mechanically-determined sparsity |15] and by 
its texture. 

Bandprofile.— Define a matrix X whose elements 
are X nm = |F„ m | 2 . The bandprofile C a (r) is obtained by 
averaging the elements X nm along the diagonals n—m = r, 
within the energy window of interest. In the same way we 
also define a median based bandprofile C s (r). See Figj3] 
for numerical results. The mean level spacing is 



A„ 



2ir/(mL x L y ) 



(8) 



Given that Ao is small compared with the energy range of 
interest, it is well known [14] that 



C a (n - m) 



2ir 



C(E n — E. m ) 



(9) 



where C(w) is the classical power spectrum, that can be 
obtained via the Fourier transform (FT) of the classical 
auto-correlation function (F(Q)F(t)). In the numerical 
analysis F(t) corresponds to a very long ergodic trajec- 
tory. It consists of impulses, namely 



F(t) = J2 2mv * 



cos 



9j)6(t-ti) 



(10) 



where 9j is the collision angle with the piston at time tj. 
By the Wicner-Khinchin theorem C(uj) oc |i^| 2 , where 
F w = FT[F{t)}. For technical details see [17] • The result 
of the calculation is displayed in Figj3] (black continuous 
line). Comparing with the quantum one observes that 
the applicability of Eq. ^ to the analysis of our billiard 
system is confirmed down to very small frequencies. 

Analytical results for C(ui) can be obtained. For large 
frequencies the power spectrum becomes flat and reaches 
the constant value [T2] 
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(11) 



For intermediate frequencies the effect of the deformation 
is mainly to ergodize the collision angle and one can ob- 
tain analytical expression (represented in Fig[3]by dashed 
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Fig. 5: The size distribution of the matrix. Histogram of 
the values of X nm for the central band of the EW1 and EW2 
matrices as defined in Fig[3] For sake of comparison we display 
results also for R = 1. 



red line) . For small frequencies the effect of the deforma- 
tion is less trivial and we find that the power spectrum is 
logarithmically divergent: 



C{u « l/i L ) » m 2 v 3 E A l n JL 
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(12) 



The divergence comes because there are vertically bounc- 
ing trajectories with very long horizontal bouncing period, 
as in the related analysis of [TB]. Disregarding the log- 
arithmic term one observe that compared with C M the 
bouncing leads to enhancement by factor 1/u, which is 
the ratio t R /i L . 

Sparsity and Texture.— For strongly chaotic sys- 
tems the elements within the band have approximately a 
Gaussian distribution. But in the WQC regime the ma- 
trix becomes sparse and textured as demonstrated in Fig[2J 
Loosely speaking, sparsity means that only a small frac- 
tion (s -C 1) of elements are largqj, while the texture refers 
to their non-random arrangement. In the WQC regime 
the size distribution of the in-band elements becomes log- 
wide (approximately log-normal) as seen in FigjSJ This is 
reflected by having 



C s (r) « C a (r) 



(13) 



as seen in Fig 01 

The sparsity and the texture of X are important for the 
analysis of the energy absorption rate [T3j as implied by 
SLRT [TQl[n]- Accordingly, we suggest to characterize 
the sparsity by a resistor network measure 



9s 



9s[X] 



((X)) s /((X)) a 



(14) 



1 A precise definition of the sparsity s can be found in Section III 
of 1191 . but it is not of much physical interest for us. Rather we char- 
acterize the sparsity by the resistor-network measure g s as defined 
below, which has direct relation to the response analysis. 
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Here ((X)) a is the algebraic average over the in-band el- 
ements of the matrix, while {{X)) s is the corresponding 
resistor network "average" that takes their connectivity 
into account. The recipe of the resistor network calcu- 
lation is detailed in the next paragraph (can be skipped 
in first reading) . For a strictly uniform matrix g s = s = 1, 
for a Gaussian matrix s = 1/3 and g s ~ 1, while for sparse 
matrix s, g s <C 1. 

The resistor network quantity {{X)) s can be regarded 
as a smart average over the elements of X, that takes 
their connectivity into account. For the purpose of its 
calculation we associate with X a matrix g whose elements 
arc 



28o(n— 7n ) 



A,, 



(n — my 



(15) 



where ^2 r o~o(r) = 1 is a weight function, whose width 
should be quantum mechanically large (i.e. ;§> 1) but 
scmiclassically small (i.e. < the bandwidth). If we take 
this weight function to be the normalized version of S(u), 
then g„ m can be interpreted as the (normalized) Fermi- 
golden rule transition rates that would be induced by a 
low-frequency driving. Optionally we can regard these 
gnm as representing connectors in a resistor network. The 
inverse resistivity of the strip can be calculated using the 
standard procedure, as in electrical engineering, and the 
result we call ((X)) s . It is useful to notice that if all the 
elements of X are identical, then ((X)) s equals the same 
number. More generally {{X)) s is smaller than the con- 
ventional algebraic average ((X)) a (calculated with the 
same weight function). In the RMT context a realistic 
estimate for ((X)) s can be obtained using a generalized 
variable-range- hopping procedure |19j . 

The WQC regime.— With the classical t L and t H , we 
can associate the energies 



A L 

Aa 



27r/t L 
2ir/U 



(16) 
(17) 



Conversely, with the mean levels spacing we can asso- 
ciate the Heisenberg time i H = 27r/Ao. Note that t H = 
(l/h) d ~ 1 t Ij where d=2. It is also possible to define the 
Ehcrnfest time t E = [log(l/ft)]£ R , which is the time re- 
quired for the instability to show up in the quantum dy- 
namics. The traditional condition for "quantum chaos" 
is t E <C £ H , but if wc neglect the log factor it is simply 
£r "C t H . This can be rewritten as A R 3> Ao, which we 
call the frequency domain version of the quantum chaos 
condition. Optionally one may write a parametric version 
of the quantum chaos condition, namely u 3> Ub, where 



Ub = ft 



[de-Broglie deformation] 



(18) 



The frequency domain version implies that it should be 
possible to resolve the zero frequency peak of C{oS) as 
in Fig.[3J while the parametric version means that a de- 
Broglie wavelength deformation of the boundary is re- 
quired to achieve "Quantum chaos" . 



Wc observe in the upper panel of FigQ]that g s is signif- 
icantly smaller than unity, even for very small values of ft 
for which u > Ub is definitely satisfied. For completeness 
we show in the lower plot additional data points in the 
regime u < Ub where this breakdown of QCC is not a big 
surprise. We conclude that QCC for u > Ub is restricted 
to C(u), and does not imply Hard quantum chaos (HQC), 
but only WQC. In the WQC regime C s (r) < C a \r) and 
consequently g s Cl, indicating sparsity. 

The emergence of WQC instead of HQC can be ex- 
plained as follows. If a wall of a billiard is deformed, the 
levels arc mixed. FOPT is valid provided \U nm \ < Ao- 
This condition determines a parametric scale u c . If the 
unperturbed billiard were chaotic, the variation required 
for level mixing would be [16] u c sa A E /(fc E L) 1 / 2 = ft 3 / 2 . 
This expression assumes that the eigenstates look like ran- 
dom waves. In the Wigner regime (u c < u < Ub) there is a 
Lorentzian mixing of the levels and accordingly, the num- 
ber of mixed levels is ~ (u/u c ) 2 . But our unperturbed 
(rectangular) billiard is not chaotic, the unperturbed levels 
of the non-deformed billiards are not like random waves. 
Therefore, the mixing of the levels is non-uniform. 

By inspection of the C/„ x „ . mxm matrix elements one 
observes that the dominant matrix elements that are re- 
sponsible for the mixing are those with large n x but small 
\n y —m y \. Accordingly, within the energy shell E n:cn ~ E, 
the levels that are mixed first are those with maximal n x , 
while those those with minimal n x are mixed last. The 
mixing threshold for the former is 



A E /(fc E L) = ft 2 



(19) 



while for the latter one finds u^ ~ ft , which is much 
larger than Ub = ft 1 . Straightforward analysis of this mix- 
ing (extending that of [12]) leads to the result 



(J 



u 2 /h 



(20) 



This is merely the ratio of the median value to the mean, 
and the proportionality to u 2 is the remnant of FOPT. 
This simple dependence is confirmed by the numerics of 
FigJU We note that the RMT perspective of [TI5] implies 
that in general this median based estimate should be cor- 
rected. Roughly the prescription is 



<! 



max{ 1 , g exp y — In b In g } 



(21) 



where b = uj c /Aq is the dimensionless bandwidth. 

In the numerics g is calculated for a bandwidth match- 
ing spectral width, i.e. the spectral support of S(u) is 
assumed to be ~ A R , implying g c ~ 0(1) and g ~ g s . In 
the quantum mechanical LRT calculation which is pre- 
sented in Figf4]by black line g c depends on ft, because A 
provides a lower cutoff on the logarithmically divergent 

<5(«). 

If the spectral support of the driving were <C A R , the 
classical correlation factor would be g c ~ 1/u, and con- 
sequently g s ~ u 3 /ft. Still, the bandwidth is the signifi- 
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cant scale in the "quantum chaos" perspective, and there- 
fore the parametric scale that signifies the WQC-HQC 



crossover is 



h 1 ' 2 



(22) 



which is larger than ut, = h. Accordingly, the WQC 
regime extends well beyond the traditional boundary of 
the Wigner regime, and in any case it is well beyond the 
FOPT border u c . 

Discussion.— In a broader perspective the term "weak 
quantum chaos" is possibly appropriate also to system 
with zero Lyapunov exponent (t K =oo), e.g. the triangular 
billiard [50], and pseudointcgrable billiards [2JJ, and to 
systems with a classical mixed phase space. But in the 
present study we wanted to consider a globally chaotic 
system, under semiclassical circumstances such that A R 
is quantum mechanically resolved and QCC is naively ex- 
pected. In this context there are of course other interest- 
ing aspects, such as bouncing related corrections to Weyl's 
law [55], and non- universal spectral statistics issues (see 
below), while our interest was with regard to the semi- 
linear response characteristics of the system. 

The spectral statistics in the WQC regime has been 
studied in [23] concerning nearly circular stadium bil- 
liard, and in [21] concerning circular billiards with a rough 
boundary. Let us remind very briefly how the WQC bor- 
der is determined in this context. It is convenient to de- 
scribe the dynamics using a Poincare map, which relates 
the angle 6 T of successive collisions (r = 1, 2, 3, • • ■) with 
the piston. One observes that due to the accumulated ef- 
fect of collisions with the deformed boundary, there is a 
slow diffusion of the angle with coefficient Dg ~ u 2 . Ac- 
cordingly the classical ergodic time is r r ~ 1/ Dg, and the 
quantum breaktime due to a dynamical localization effect 
is Th ~ De/ti 2 . The border of the WQC regime is defined 
by the condition Th < t t leading to Eq. (f2"2"|) . However we 
would not like to over-emphasize this consistency because 
it is not a-priori clear that spectral-statistics and sparsity 
related characteristics always coincide. 

Practical implications.— Coming back to the "con- 
flicting expectations" issue, with regard to the value of the 
absorption coefficient and its dependence on the deforma- 
tion u, we now can see how they reconcile. First of all it 
should be clear that if there were no classical correlations 
between bounces, then C(ui) would be flat, equals to the 
value Coo of Eq. ([TT)) . leading to the wall formula Eq. flBJ) 
for G. The effect of bouncing is to enhance C(ui <C 1/*l) 
as implied by Eq. ([T2"]) . Depending on whether the spectral 
support of the driving is oj c <C A R or u c ~ A R we observe 
or do not observe al/u enhancement. This holds classi- 
cally and also in the quantum LRT calculation (provided 
uj c > A ) due to QCC. 

However, the SLRT calculation, unlike the LRT calcu- 
lation, cares about the median and not about the mean. 
Therefore, for a weakly chaotic system, it give a much 



smaller result for G. If the mixing of the levels were uni- 
form we would expect a crossover from g s oc u 2 (FOPT) to 
g s ex 1/u 2 (Wigner), as in the theory of disordered conduc- 
tors. But the mixing of levels in a weakly chaotic system, 
unlike in a weakly disordered system, is not uniform, and 
therefore the g s oc u 2 persists within a very large range 
u c < u < u s , to which we refer as the WQC regime. 

Experimental feasibility.— Having a better under- 
standing of the WQC regime we are now able to revise the 
suggested experiment in [TJ]. Let us consider S5 Rb atoms 
that are laser cooled to low temperature T w 0.1/iK , 
such that the de-Broglie wavelength is A E = 1 fim. The 
atoms are trapped in an optical billiard of linear size of 
L = 10 fim, and accordingly the dimensionless Planck con- 
stant is h = 0.1. This leads to A L /A = 30. Note that 
A L = 220 Hz, and A = 7.5 Hz. 

Assuming 10% deformation the dimensionless band- 
width can be tuned as b = (A R /A ) ~ 10. By modulating 
the laser intensity, one of the billiard walls can be noisily 
vibrated. We assume that the driving is band-matching, 
i.e. uj c ~ A n . These are roughly the same parameters as 
in our numerical analysis. The prediction for the SLRT 
suppression factor is g s ~ 0.1. 

In order to witness the SLRT anomaly the RMS ampli- 
tude of the vibrations (e) should be large enough, as to 
have a measurable heating effect. Assuming that it is pos- 
sible to hold the atoms for a duration of ~ 1000 bounces 
the condition can be written in a dimensionless form as 
G £ 2 /(TA L ) > 10- 3 , or roughly as (e/L) 2 > 10~ 3 . 

On the other hand e should be small enough, such 
that the FGR condition is not violated. It is straight- 
forward to show that the FGR condition can be written 
in a dimensionless form as TGqE 2 / A 3 , < b 3 , or roughly 
as (e/L) < (1/6). Accordingly there is a range where 
both conditions are satisfied, and there the SLRT anomaly 
should be observed, provided environmental relaxation ef- 
fects can be neglected. 

Comments.— It is important to realize that we are 
studying in this work a driven chaotic system, and not a 
driven integrablc system. Remarkable examples for driven 
integrablc systems are the kicked rotator [25] and the vi- 
brating elliptical billiard [2E] . In the absence of driving 
such systems are integrable, while in the presence of driv- 
ing a mixed phase space emerges. This is not what we call 
here weak chaos. 

The low frequency driving that we assume is stochastic, 
rather than periodic. This looks to us realistic, reflect- 
ing the physics of cold atoms that are trapped in optical 
billiards with vibrating walls. It is also theoretically con- 
venient, because we can use the FGR picture. If one is 
interested in periodic driving of strictly isolated system, 
then there are additional important questions with regard 
to dynamical localization [57], that can be handled e.g. 
within the framework of the Floquet theory approach. 
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Summary.— The discovery of "anomalies", i.e. ma- 
jor deviations from QCC in circumstances where QCC 
is expected by common wisdom, is a major challenge in 
quantum-mechanics studies. For example: Anderson's Lo- 
calization (wavefucntions were commonly expected to be 
extended); Heller's scars (wavefucntions were commonly 
expected to look like random waves) . Here we highlighted 
an anomaly in the theory of response: the rate of heating 
is unexpectedly suppressed for a quantized chaotic system. 

Our analysis has been based on SLRT. This theory ap- 
plies to circumstances in which the environmental relax- 
ation is weak compared with the /(i)-induced transitions. 
In such circumstances the connectivity of the transitions 
from level to level is important, and the LRT result should 
be multiplied by g s . 

We have highlighted that there is a distinct WQC 
regime, where semiclassics and Wigner-type mixing co- 
exist. This is the regime where an LRT to SLRT crossover 
is expected as the intensity of the driving is increased. 

Acknowledgements.— We thank Nir Davidson (Weiz- 
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Binational Science Foundation (BSF). 



REFERENCES 



[I] D.H.E. Gross, Nucl. Phys. A 240, 472 (1975). J. Blocki, Y. 
Boneh, J.R. Nix, J. Randrup, M. Robel, A.J. Sierk, W.J. 
Swiatecki, Ann. Phys. 113, 330 (1978). S.E. Koonin, R.L. 
Hatch, J. Randrup, Nuc. Phys. A 283, 87 (1977). 

[2] N. Friedman, A. Kaplan, D. Carasso, N. Davidson, Phys. 

Rev. Lett. 86, 1518 (2001). 
[3] A. Stotland, R. Budoyo, T. Peer, T. Kottos, D. Cohen, J. 

Phys. A 41, 26200RFTC) (2008). 
[4] E. Ott, Phys. Rev. Lett. 42, 1628 (1979). R. Brown, E. Ott, 

C. Grebogi, Phys. Rev. Lett. 59, 1173 (1987). R. Brown, 

E. Ott, C. Grebogi, J. Stat. Phys. 49, 511 (1987). 
[5] C. Jarzynski, Phys. Rev. E 48, 4340 (1993). C. Jarzynski, 

Phys. Rev. Lett. 74, 2937 (1995). 
[6] M. Wilkinson, J. Phys. A 21, 4021 (1988). M. Wilkinson, 

E.J. Austin, J. Phys. A 28, 2277 (1995). 
[7] J.M. Robbins, M.V. Berry, J. Phys. A 25 L961 (1992). 
[8] D. Cohen, Phys. Rev. Lett. 82, 4951 (1999). D. Cohen, 

Annals of Physics 283, 175 (2000). D. Cohen, T. Kottos, 

Phys. Rev. Lett. 85, 4839 (2000). 
[9] D.M. Basko, M.A. Skvortsov, V.E. Kravtsov, Phys. Rev. 

Lett. 90, 096801 (2003). A. Silva, V.E. Kravtsov, Phys. 

Rev. B 76, 165303 (2007). 
[10] D. Cohen, T. Kottos, H. Schanz, J. Phys. A 39, 11755 

(2006). 

[II] M. Wilkinson, B. Mehlig, D. Cohen, Europhys. Lett. 75, 
709 (2006). 

[12] A. Stotland, D. Cohen, N. Davidson, Europhys. Lett. 86, 

10004 (2009). 
[13] R. Ram-Mohan, Finite Element and Boundary Element 

Applications in Quantum Mechanics (Oxford University 

Press, Oxford, UK, 2002). 



[14] M. Feingold, A. Peres, Phys. Rev. A 34 591, (1986). M. 

Feingold, D. Leitner, M. Wilkinson, Phys. Rev. Lett. 66, 

986 (1991). 
[15] E.J. Austin, M. Wilkinson, Europhys. Lett. 20, 589 (1992). 

T. Prosen, M. Robnik, J. Phys. A 26, 1105 (1993). Y. 

Alhassid, R.D. Levine, Phys. Rev. Lett. 57, 2879 (1986). 

Y.V. Fyodorov, O.A. Chubykalo, F.M. Izrailev, G. Casati, 

Phys. Rev. Lett. 76, 1603 (1996). 
[16] D. Cohen, A. Barnett, E.J. Heller, Phys. Rev. £63, 46207 

(2001). 
[17] A. Barnett, D. Cohen, E.J. Heller, Phys. Rev. Lett. 85, 

1412 (2000); J. Phys. A 34, 413 (2001). 
[18] B.J. Alder, T.E. Wainwright, Phys. Rev. A 1, 18 (1970). 

F. Vivaldi, G. Casati, I. Guarneri, Phys. Rev. Lett. 51, 727 

(1983). 
[19] A. Stotland, T. Kottos, D. Cohen, Phys. Rev. B 81, 

115464 (2010). 
[20] G. Casati and T. Prosen, Phys. Rev. Lett. 85, 4261 (2000) 

M. Degli Esposti, S. OKeefe and B. Winn, Nonlinearity 18, 

1073 (2005). 
[21] E.B. Bogomolny, U. Gerland, C. Schmit, Phys. Rev. E 59, 

R1315 (1999). 
[22] A. Backer, R. Schubert, P. Stifter, J. Phys. A 30 6783 

(1997). 
[23] F. Borgonovi, G. Casati and B. Li, Phys. Rev. Lett. 77, 

4744 (1996). 
[24] K.M. Frahm and D.L. Shepelyansky, Phys. Rev. Lett. 78, 

1440 (1997). 
[25] B.V.Chirikov, Phys. Rep. 52, 263 (1979). 
[26] F. Lenz, F.K. Diakonos, P. Schmelcher, Phys. Rev. Lett. 

100, 014103 (2008); Europhys. Lett. 79, 2002 (2007). 
[27] T. Prosen, D.L. Shepelyansky, Eur. Phys. J. B 46, 515 

(2005). S. Fishman, D.R. Grempel and R.E. Prange, Phys. 

Rev. Lett. 49, 509 (1982). 



p-7 



